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THE MINIMUM FREE-STREAM WIND SPEED FOR INITIATING MOTION 
OF SURFACE MATERIAL OE MARS 

By G. P. Wood, W. R. Weaver, and R. M. Henry 
Langley Research Center 

SUMMARY 

Estimates of the minimum free-stream wind speed that is required for 
initiating the notion of surficial material on Mars have ranged from 30 to 
about 200 Refers per second. Thus the best value for this quantity is not well 
established. Graphical comparison of such of the pertinent data taken in the 
laboratory and in the field cn Earth provides a minima value for the "Bagnold 
coefficient" of 0.08 and this in turn provides a minimum value for the thresh- 
old friction velocity of 1.3 m/s for initiating motion of particulate matter on 
Mars at low elevations where the pressure is 7 mb. The most appropriate value 
of the ratio of friction velocity to free-stream velocity for putative unstable 
conditions appears to be 0.026. Thus the minimum free-stream wind apeed for 
initiating motion is obtained as SO m/s. If the surface material on Mars, 
however, is less cohesive than that on Earth, the minimum value may be smaller. 

INTRODUCTION 

The minimum free-stream wind speed required to initiate the motion of sand 
or dust on the surface of Mars has been estimated in a number of papers, with 
results that have been interpreted as ranging from 30 to perhaps as high as 
200 m/s. Sagan and Pollock (ref. 1) are sometimes quoted as deriving the mini- 
mum required free-stream speed as being 80 to 200 m/s for a surface pressure of 
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S a to. This result can easily be obtained from their paper, but it would be 
fairer to quote a different result fro* that paper. They state that, for a 
pressure of 10 ab, if the atmosphere is neutrally stable the required speed is 
in excess of 83 */s but, since the ataospher ; is highly unstable during the day, 
■ilder winds than that would be required. But they add that the expected 
uncertainties are factors of two or three. Golitsyn (ref. 2) estimates the 
ainiaim required speed to be about 80 m/s. Hess (ref. 3) obtains, for a surf- 
ace pressure of 8 ab, a free- stream speed of 60 a/s. Gierasch and Goody (ref. 

4) derive the surpiisingly low value of 30 a/s whereas Iversen et al. (ref. 5) 
have the surprisingly high estimate of 120 a/s. 

Kith the infoxaation that is available, it appears to be possible to make 
a better estiaate of the value of th* ainiaua required free-streaa wind speed 
than has soaetiaes been done. It is the purpose of this paper to sake this 
estiaate. 


SYMBOLS 

A coefficient in Bagnold's equation for u* 

t 

A,B quantities in equation for ratio of velocities 

C D proportionality constant 

D particle diaaeter 

f Coriolis parameter, 2u sin $ 

g acceleration due to gravity 

Re* particle friction Reynolds number 

Ro surface Rossby number, U/fz Q 

U free-streaa wind speed 

1/2 

u* friction velocity, (r/p) 
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threshold value of friction velocity 

roughness length 

von Karman constant 

coefficient of viscosity 

fluid density 

particle density 

shear stress 

latitude 

angular speed of rotation of planet 


THE THRESHOLD VALUE OF FRICTION VELOCITY 


The stress T exerted on a surface by a wind is proportional to the den- 
sity p and to the square of the speed of the wind at the top of the boundary 
layer. 


x = C D pU 2 . 

If one defines a speed u„, which is generally called the "friction velocity," 
as Cp*^ 2 U, then 


2 , 
u * = T/p 

It has been shown (e.g., by Bagnold, ref. 6) that the threshold value of the 
friction velocity, which is the value at which the wind initiates motion of 
grains of diameter D, can be calculated hy the equation 


u* 


t 





where, for cohesionless grains, the coefficient A is a function of the 
"particle friction Reynolds number" 
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pu.D 

Re * * “ 


( 2 ) 


For a turbulent boundary layer (above the laminar sublayer, if a sublayer 
exists) A is a veiy weak function of Re*. For the laminar case, A is a 
very rapidly varying function of Re*. Most of the data on threshold condi- 
tions have been presented in the form of plots of threshold value of friction 

velocity u* versus effective diameter of particles. In the extension of 
t 

terrestrial data to Mars, it has been assumed by others (e.g., Ryan, ref 7, 

and Gifford, ref. 8) that A is a universal function of Re*. To facilitate 

this extension, we have plotted in figure 1 the results of six papers (refs. 5, 

9 to 13) in the form u* as a function of Re* for terrestrial conditions 

t 

(density, gravity, viscosity), he have made the same assumption as others and 

have converted these results to the form of the Bagnold coefficient A as a 

function of Re* in figure 2, for application to Mars. 

The smallest values of u* shown on figure 1 are 0.14 m/s for Bagnold 

t 

(ref. 9) and 0.12 and 0.10 for Chepil (refs. 11 and 12). For the first of 

these the nominal diameter represents a range of diameter for which the ratio 

of maximum to minimum apparently was 1.8. The second was for soil with 15% 

clods of up to 25-mm size, and the third was for a range of sizes from fine 

dust up to the nominal size. Chepil in reference 12 implies that for a mixture 

of grains "such as a commonly occurring dune material," the coefficient A is 

reduced by about 15% at the average equivalent diameter of the mixture. Thus, 

for a mixture that is composed of a narrow range of sizes, we would use 0.14 

for u* . For a mixture composed of a wide range of sizes, such as dune sand, 
t 

we feel justified in using a value as small as 0.12, and this is the value that 
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we have chosen. The corresponding value of Re* is seen to he 0.9. From 
figure 2 the corresponding value of A is 0.08. 

With values for A and Re*, the two equations (1) and (2) can he solved 
for the minimum value of u* and the corresponding value of D. We use a 
particle density a of 2,65 x IQ 3 kg/m 3 and a coefficient of viscosity y of 
1.1 x io S kg/m sec. We use atmospheric densities corresponding to two eleva- 
tions on Mars: altitude zero with respect to the mean surface level in the 

-2 3 

equatorial zone, for which the density is 1.2 x 10 kg/m and the pressure is 

_2 

4.8 mb, and altitude 4 km below that level, for which the density is 1.7 x IQ 

3 

kg/m and the pressure is about 7 mb. This latter altitude is approximately 

the altitude of the proposed Viking landing sites in Chryse and Tritonis Lacus 

(formerly Nodus Laocoontis). Use of the altitude of these two landing sites is 

particularly appropriate. They may be the locales of relatively frequent events 

or changes in appearance that can be attributed to wind-blown sand or dust. 

Antoniadi's book (ref. 14) states that yellow clouds very frequently obscure 

Chryse. Slipher’s book (ref. 15) states that ’The large bright region between 

Elysium and Syrtis Major has revealed the greatest array of striking changes 

(between 1907 and 1962J of any part of the desert areas of Mars Perhaps 

the most extraordinary change in this region occurred around Nodus Laocoontis...." 

Simultaneous solution of equations (1) and (2) with A = 0.08 and 

Re* =0.9 provides the results that for altitude zero the minimum of the 

threshold value of the friction velocity u* is 1.6 m/s and the corresponding 

t 

value of diameter D is 510 yin; for the lower altitude the minimum value of 

u* is 1.3 m/s and the corresponding value of D is 450 ym. 
t 
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POSSIBILITY OF REDUCED COHESION 


It must be recognized that the process of deriving a minimum value of u* 

t 

for Mars as above on the basis of data taken on Earth may have omitted an 
effect that could conceivably alter to a considerable extent the values of u* 

t 

The Bagnold relation, equation (1), was derived without taking account of possi 
ble cohesion between adjacent grains of soil. (The word "cohesion" is used 
herein to mean either or both cohesion and adhesion.) Almost all experiments 
were done without regard to the constitution of the ambient gas or its moisture 
content. It is conceivable that the CO^ atmosphere on Mars may affect the 
cohesion of the particles of the soil, and it is much more likely that the 
greatly reduced water vapor content of the Martian atmosphere reduces signifi- 
cantly any forces of cohesion among particles because of reduced quantities of 
adsorbed water. The only experiments in which the relative humidity of the air 
was controlled seem to be those reported in Addendum III of reference 16 hy 
Belly. These experiments appear to show that the amount of water vapoi in the 

air can have a large effect on u* even for grains as large as those in a 

t 

mixture having a mean diameter of 440 pm. Thus, we feel constained to point 

out that our minimum value of u* for Mars may in fact be an overestimate of 

t 

the minimum, and that the true minimum, if one exists, may be considerably 
smaller. Furthermore, if cohesion does play an important role, then it would 
appear that the general ly -made assumption that the Bagnold coefficient A is a 
function only of Re* is not supportable. 
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THE EFFECT OF FIXED ROUGHNESS ELEMENTS 


The presence of immobile, or fixed, roughness elements (usually called non- 

erodible roughness elements) can have a rather large effect on the value of 

u* . There is not, however, very much quantitative information on the effect, 
t 

Chepil has stated (ref. 17) that "where the roughness elements .. .are non- 

erodible the value of coefficient A is increased considerably." Iversen et 

al. have shown (ref. 5) that the value of the coefficient A can be tripled by 
the presence of fixed roughness elements. Chepil however shows (ref. 11) no or 
very little effect of fixed roughness near the minimum part of the curve of 

u* vs Sd. In the present calculations, we assume either an absence of fixed 

t 

roughness elements or that they do not have significant effect on the minimum 

value of u„ . 

t 

MINIMUM THRESHOLD VALUE OF THE FREE-STRF.AM SPEED 

We obtain the estimate of the free-stream speed appropriate for our value 

of u„ by determining the value of the ratio u* /l). This ratio can be con- 
t t 

sidered to be a function of the Rossby number and of the degree of stability or 
instability. The ratio can be taken to be (Shir, ref. 18) 



where Awl. 3, B « 2.8 for the neutrally stable case. This equation appears 
to be valid also for the unstable case, but the values of A and B as func- 
tions of the degree of stability are not known. 
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We use 0.40 for k, which is the value that has been in use for about 30 
years (ref. 19). Reference 19, for example, states, however, that for Rossby 
numbers corresponding to smooth terrain 0.35 is probably a good choice for the 
value of k. But in the absence of much supporting evidence for this value, we 
have used 0.40. If the correct value is less, than U would be correspondingly 
larger. 

9 

We estimate the surface Rossby number U/fz^ to be about 2 * 10 . From 

u* 

the preceding equation or from the plot of -g- vs Ro, figure 3 of Shir (ref. 18), 
we find u*/U = 0.024. If we use figure 4 of Csanady (ref. 20) we obtain 0.021. 
If we use figure 3 of Monin (ref. 21) we have 0.025. All of these results are 
for the neutrally stable boundary layer. Csanady's figure 4 does show also 
results for the stable and the unstable cases, but we believe that these are 
correct only qualitatively, not quantitatively, because the value of the constant 
in the first of Csanady's equation (19a) applies only to the neutrally-stable 
case (Hess and Clarke, ref. 22). For the unstable case the value of the constant 
should be less. Csanady's figure 4 does show that for the unstable boundary 
layer the value of u*/U is greater than for the neutrally-stable case, unless 
Ro is very large. Monin states that measurements on Earth show that as strati- 
fication increases from stable to lapse, u„/U increases. Pasquill (ref. 23) 
quotes results of Clarke (ref. 24) as showing that, on the average, for unstable 
conditions the value of u*/U is 1.6 times its value for stable conditions. 

These quantitative results, however, must be for smaller values of Ro than are 
appropriate for Mars. Therefore, it does not appear to be either unreasonable 
or unconservative to use a value of at least 0.026 for the ratio for the unstable 
case on Mars. Thus, for the higher altitude, for which u* - 1.6 m/s, the free- 
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stream wind speed U = 62 m/s. For the lower altitude, for which u* ■ 1.3, 

U » 50 m/s. 

We have thus arrived at a minimum value for the required free-stream wind 
speed for initiating the motion of surface material at the proposed landing 
sites for the 1975 Viking mission to Mars. Although we used rather different 
values than Hess did (ref. 3) for most of the quantities that enter the calcu- 
lation, we obtain a result fairly close to his. This result should allay some 
of the concern that is generated by estimates that much stronger winds are 
required. 
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Figure 1. - Threshold friction velocity for Earth's atmosphere as calculated 
from data of various experimenters. 




Figure 2. - Coefficient A as calculated from data of various experimenters 
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